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C^ , We consider a two dimensional semiconductor with carriers subject to spin- 

orbit interactions and scattered by randomly distributed magnetic impurities. 
We solve the time-dependent Schrodinger equation to investigate the relation- 
f^J ' ship between the geometrical properties of the wavefunction and the system's 

^H i spin dependent transport properties. Even in the absence of localized states, 

Q i interference effects modify the carrier diffusion, as revealed by the appearance 

r\ ■ of power laws dependent on the impurity concentration. For stronger disorder, 

we find a localization transition characterized by a fractal wavefunction and 
enhanced spin transport. 

^ . Keywords: Condensed matter; diluted magnetic semiconductors; spintronics; 

00 ' localization. 

1. Introduction 

CO ! 

f— s , An electron gas confined in a semiconductor heterostructure can display, 

{N| ■ under the action of an external electric field, spin dependent transport 

driven by the spin-orbit interaction. 1 Although this intrinsic spin Hall effect 
is suppressed by disorder, 2 magnetic impurities can restore it. 3 

ISJ{ ■ Indeed, in the absence of disorder, an electric field in the x-direction, 

^ \ induces an out of plane spin polarization following the ^/-direction. The spin 

current j* is related to the applied electric field E x by the spin conductivity 
jy = cj s hE Xi whose value cj s h = — e/87r (e is the electric charge), turns out 
to be independent of the spin-orbit coupling constant A. 1 Non-magnetic 
disorder preserves the time reversal symmetry of the microscopic Hamil- 
tonian, and allows the setting of a macroscopic stationary state. Due to 
the particular form of the Rashba term, linear in the momentum, the spin 
current is proportional to the time derivative of the in plane spin; as this 
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derivative must vanish at equilibrium, the spin Hall effect is canceled. 4 In 
the presence of magnetic disorder, the system looses its microscopic time 
reversal invariance, allowing the emergence of the intrinsic spin Hall effect, 
even if this invariance is statistically restored. 3 ' 5 

Our purpose here is to investigate the behavior of the spinor wavefunc- 
tion when the impurity concentration increases and to relate its dynamical 
and geometrical properties to the transport properties of the system. We 
present the model in Sec. 2, followed be an account of the numerical results 
(Sec. 3) and the conclusions (Sec. 4). 

2. Model 

The Hamiltonian of the two-dimensional system writes 

p 2 X 

H = — --a-(zxp)- J s a- n(x) , (1) 

where the second term represents the Rashba spin-orbit coupling, with pa- 
rameter A, and the last term describes the exchange interaction, with ex- 
change energy J s and magnetic impurities z, 

n(x) = a 2 2_. Tii5{x — xi) , (2) 

oriented in the direction n^, and distributed randomly; here X is the set of 
impurity sites in the plane x = (x, y) (a is a length constant); cr is the vector 
of Pauli matrices, m the carrier's effective mass, and p = (p x ,P y ) their 
momentum. We consider the magnetic impurities in a paramagnetic state, 
therefore the average of the last term in (1), over the magnetic moment 
orientations n^, vanishes. 

The dynamics of the quantum state \ijj) is governed by the Schrodinger 
equation, 

i^- t m)) = H\m), (3) 

where we chose units such that h = a = m = 1. We numerically integrate 
(3) using a Chebyshev spectral method, 6 to obtain the time evolution of an 
initial spin-up gaussian wave packet. We measure the width, 

wl(t)= I <te|*| 2 hfc(*,*)| 2 , (4) 

where ip(x,i) = (x\ip(t)) = (^ + ,^_) T is the spinor whose components are 
the spin up and down amplitudes, and L 2 is the system's surface. The 
initial orientation of the spin is not relevant, scattering off impurities is 



February 26, 2013 19:43 



WSPC - Proceedings Trim Size: 9in x 6in cctll 




Fig. 1. Time evolution of the absolute value of the down-spin wavefunction ip—, for 
t -1 = 0.004, 0.1. Extended state (top); nearly localized state (bottom). The four panels 
of each row are for times 20, 60, 160, and 400, respectively. 



statistically isotropic, so that asymptotically both spin states display the 
same behavior. The spin current in the (arbitrary) direction y, is given by 



1 



f y (t) = j(m\Wz,vy}\m) 



i 



-Im 



fdx(^(t)d v Mt) - r-{t)d v ^-{t)) , 

(5) 



using the spinor in the momentum representation; here v y is the y com- 
ponent of the velocity operator, and the brackets {•,•} are for the anti- 
commutator. The physical dimension is [j*] = [hp/mL 2 ], the same as a 
force: charge (coming from the spin conductivity), multiplied by an elec- 
tric field (driving force). In the absence of an external electric field, the 
transverse spin current should vanish in a stationary state, however in the 
present initial value problem, the |^(0)) is not an eigenstate of i7, leading to 
an unbalance of spin- up and spin-down probability currents, which results 
in a (at least transitory) net spin current. 

To characterize the geometry of the wavefunction we measure the inverse 
participation ratio P(t)J and the fractal dimension D(t). 8 It is convenient 
to partition the surface L 2 into Ni cells of size I 2 and to define the measure, 9 



m k 



(i,t)= J2 hM*)l 2 > fc = i,-..,w, 



(6) 



i=lek 



Q2 k rrik = 1), from which we write the dynamical ^-inverse participation 
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Fig. 2. Spin current as a function of time, corresponding to r 1 = 0.004 (left) and 
t -1 =0.1 (right), the two cases of Fig. 1. 



ratio, 



N L 



P q (l,t) = Y t m k (l,ty, 



(7) 



fe=i 



1 log P q (l,t) 



(8) 



and q- fractal dimension, 10 

DM) = lim- , , 

gW z^og-1 log/ 

We use in the following the special case with q = 2, D = D2 and P = P2. 
The wavefunction fractal dimension D = D2, determines the asymptotic 
behavior of the correlation functions, or return probability, ~ f-D 2 /2 j n 
the critical state. 11 The inverse participation ratio tends to zero when the 
wavefunction is extended, and remains finite (at long times) when the quan- 
tum state is localized; concomitantly, the fractal dimension drops from a 
value of 2 (the entire surface) in the extended regime, to a smaller value in 
the localized regime. These quantities are well suited for characterizing the 
Anderson transition as a function of the disorder. 12 The relevant param- 
eter to characterize the disorder strength is the scattering time r, which 
is proportional to the exchange energy and the impurities concentration c, 



3. Results 

We analyze the temporal evolution of an initially normalized Gaussian wave 
packet, concentrated around the origin, 

1 / Jaf\ 



M^,o) = 



(27rd 2 y/- 



exp 



4d 2 /' 



(9) 
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Fig. 3. Spreading of the wavefunction t/l_ for r~ 1 = 0.004 (left) and r~ 1 = 0.1 (right), 
as in the previous figures. They follow a power law with exponent (3 (slope of the solid 
line). 



and ip-(x : 0) = 0, where the size is typically d = 12a, much smaller than the 
system size (^ 512a). We vary the value of the disorder strength, essentially 
by changing the exchange coupling constant. In a clean system, which can 
be taken as a reference, the width of the wave packet spreads ballistically, 
w ~ t, its correlation dimension D = 2, and its inverse participation ratio 
tends to zero. 

In Fig. 1 we compare the time evolution of the Gaussian wave packet 
for two values of the disorder strength: r _1 = 0.1 x (0.2) 2 , 0.1 x (1.0) 2 (the 
impurity concentration is c = 0.1). We note, on the one hand, that in the 
weak disorder case (top) the continuous spreading of the wavefunction is 
accompanied by a significant amount of interference; on the other hand, for 
stronger disorder (bottom), the spatial distribution of the wavefunction is 
more complex, and the spreading slower. Associated with these two regimes, 
we observe that the spin current (for one realization of disorder) is greatly 
enhanced in the strong regime, as shown in Fig. 2. The order of magnitude 
of the amplification factor is the ratio of the respective scattering strengths. 

This somewhat paradoxical result means that, increasing the disorder 
and consequently decreasing the electron mobility, the spin current (fluc- 
tuations) intensifies. Indeed, the spin current is not directly related to the 
(ballistic) mobility of the carriers, but to their drift "velocity" (spin tilting) 
in an effective, time dependent, magnetic field created by the cyclotron 
spin-orbit motion. We remind that in our case, the spin current result from 
the time dependence of the quantum state, together with the vanishing of 
the probability current, and hence the charge current (the spin current is 
invariant under time reversal, at variance with the charge current). 
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Fig. 4. Fractal dimension for r 1 = 0.004 (left) and r 1 =0.1 (right), as in the 
previous figures. The strong disorder case is near the localization transition. 



These qualitative features are more quantitatively characterized by the 
measure of the wavefunction width (Fig. 3), the fractal dimension (Fig. 4), 
and the inverse participation ratio (Fig. 5). We find the the asymptotic 
behavior of the width is well described by a power law w(i) ~ t& . However, 
in contrast to the usual ballistic /3 = 1 or diffusive behavior /3 = 0.5, the ex- 
ponent explicitly depends on the disorder strength, as seen in Fig. 3, where 
f3 = 0.84 and 0.46. This is confirmed by the measure of the fractal dimension 
of Fig. 4, that initially increases slowly with time, but remains well below 2 
(the value of the extended state, filling the surface) in the strong disorder 
case. The slower spreading rate can be attributed to a sticking effect of the 
wavefunction around the impurities, leading simultaneously to a complex 
spatial distribution (as reflected by the smaller correlation dimension). 

We may suppose that this later case should be near a localization tran- 
sition. In order to verify this prediction, we compare the behavior of the in- 
verse participation ratio for three different values of the scattering strength, 
the two used so far, and a very strong one, r _1 = 2.5, for which we must 
certainly be in the localized state. We plot P(t) in Fig. 5, for the three 
cases; we note that already the intermediate case, shows a saturation to 
a finite value, indicating that the wavefunction is actually localized. The 
fractal dimension in this case tends to the stationary value of D = 1.33 (c.f. 
Fig. 4). In the very strong disorder case, the fractal dimension saturates at 
D = 1.0. 

Therefore, around the value r _1 ~ 0.1 a localization transition takes 
place, characterized by are a drastic reduction in the wave packet spread- 
ing rate, an important enhancement of spin current fluctuations, and a 
correlation dimension that approaches one. 
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Fig. 5. Inverse participation ratio, for three scattering times, r _1 = 0.004 (left), 0.1 
(center), and 2.5. Note the different time ranges. The intermediate case is already in the 
localization regime. 

4. Conclusion 

We demonstrated that the transport properties, diffusion of the wavefunc- 
tion, spin current, in a two-dimensional semiconductor heterostructure are 
related to the localization of the quantum state. The interplay of spin-orbit 
coupling and magnetic disorder leads to atypical behavior of the wavefunc- 
tion spreading, with a power law that depends on the scattering time in the 
extended regime. These dependence can be attributed to the importance of 
interference processes, that partially localize the wavefunction around the 
impurities. 

For strong disorder, the system undergoes an Anderson transition. The 
localization of the quantum states is accompanied by the appearance of 
strong fluctuations in the spin current. In addition, the wavefunction devel- 
ops a fractal structure with a correlation dimension well below the normal 
value, corresponding to the system's space dimension, reflecting a highly 
intermittent probability distribution. These spatial fluctuations allow the 
existence of large, local spin up and down fluctuations, contributing to spin 
transport. 

It would be worth studying the behavior of the local density of states 
and its relation with the spin Hall conductivity. This investigation is in 
progress. 
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